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Abstract. We propose a generalized accretion disk viscosity 
prescription based on hydrodynamically driven turbulence at 
the critical effective Reynolds number. This approach is con- 
sistent with recent re-analysis by ^chard & Zahn (1999) of ex- 
perimental results on turbulent Couette-Taylor flows. This new 
/3-viscosity formulation is applied to both selfgravitating and 
non-selfgravitating disks and is shown to yield the standard a- 
disk prescription in the case of shock dissipation limited, non- 
selfgravitating disks. A specific case of fully selfgravitating 
/3-disks is analyzed. We suggest that such disks may explain 
the observed spectra of protoplanetary disks and yield a natural 
explanation for the radial motions inferred from the observed 
metallicity gradients in disk galaxies. The /3-mechanism may 
also account for the rapid mass transport required to power ul- 
tra luminous infrared galaxies. 

Key words: Accretion, accretion disks - Hydrodynamics - 
Turbulence - Stars: pre-main sequence - Galaxy: evolution - 
galaxies: evolution 



1. Introduction 

One of the major shortcomings of the current theoretical de- 
scriptions of accretion disks is lack of detailed knowledge 
about the underlying physics of viscosity in the disk. This prob- 
lem is significant because almost all detailed modelling of the 
structure and evolution of accretion disks depends on the value 
of the viscosity and its dependence on the physical parame- 
ters. There is general agreement that molecular viscosity i^moi 
is totally inadequate and that some kind of turbulent viscosity 
is required. Moreover, the Reynolds number in the disk flow is 
extremely high in any astrophysical context and this in itself is 
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likely to lead to strong turbulence regardless of the details of 
the actual instability involved. 

However, there is far less certainty about how to prescribe 
such a turbulent viscosity in the absence of a proper physical 
theory of turbulence. Most inv e stigat ors a dopt the so-called a 



ansatz introduced by |Shakura| ( |1972| ) and [Shakura & Sunyaev 



( |1973| ) that gives the viscosity {v) as the product of the pressure 
scale height in the disk {h), the velocity of sound (cg), and a 
parameter a that contains all the unknown physics: 



= ahc 



(1) 



One interprets this as some kind of isotropic turbulent viscosity 
V = vt = ItVt where It is an (a priori unknown) length scale 
and vt an (a priori unknown) characteristic velocity of the tur- 
bulence. One may then write a = (vt/cs) • {h/h). On general 
physical grounds neither term in parentheses can exceed unity 
so that a < 1. If initially vt > Cg, shock waves would result 
in strong damping and hence a return to a subsonic turbulent 
velocity. The condition It > h would require anisotropic turbu- 
lence since the vertical length scales are limited by the disk's 
thickness, which is comparable to h. 

While it is always, in a trivial way, possible to calculate a 
value a, a parameterization of this sort for f is only useful if 
the proportionality parameter, a, is (approximately) constant. 
One can expect this to happen only if the scaling quantities 
are chosen in a physically appropriate manner. Models for the 
structure and evolution of accretion disks in close binary sys- 
tems (e.g., dwarf novae and symbiotic stars) show that Shakura 
& Sunyaev's parameterization with a constant a leads to re- 
sults that reproduce the overall observed behaviour of the disks 
quite well. Time dependent model calculations of the outbursts 
of dwarf novae (e.g., ^eyer & Meyer-Hofmeister 1984) and 
X-ray transients (e.g., Zannizzo 1996) demonstrate that, over 



a wide range of physical states of a disk in different phases of 
its evolution, the derived values of the viscosity parameter a 
do not vary by more than approximately an order of magnitude 
and are not too different from unity. As a result of this suc- 
cess, the a-ansatz is now used in essentially all accretion disk 
models. 

It is noteworthy that, despite this success, the a-ansatz re- 
tains no information about the mechanism generating the turbu- 
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lence but only about physical limits to its efficiency in a disk. In 
fact we would expect any high Reynolds number astrophysical 
shear flow to exhibit some kind of turbulent viscosity regardless 
of whether or not it happens to be in a disk. We therefore con- 
clude that a more general prescription underlies the a-ansatz 
for accretion disks. 



In recent years, Balbus & Hawley (1991) and their collabo- 
rators (e.g., Hawley, Gammie & Balbus 1995) have shown that 
for non-selfgravitating magnetic accretion disks, an instability 
exists that can give rise to turbulence with the required formal 
dependence and — if only marginally — the required amount. 
We also note that in substantial regions of proto-stellar and 
proto-planetary disks, the charge density is unlikely to be high 
enough to sustain a significant magnetically mediated viscos- 
ity, although this phenomenon may be relevant elsewhere. 

Whether a purely hydrodynamic turbulence can sustain the 
viscosity in the angular momentum profile of an accretion disk 
and can result in an angular momentum transport towards re- 
gions with larger specific angular momentum is still a matter 



of debate. Balbus, Hawley & Stone (1996), for instance, argue 



against it, based on numerical experiments, albeit for a rather 
low effect ive Reynolds number. [Dubrull^ (|1992|) and K ato & 
Yoshizawa ( |1997| ), among others, argue in favor of it, mainly 
based on analytical considerations. Experiments dating back to 
the 19 30s on the Couette-Taylor fl ow between co-axial rotating 
fluids ( Wendt 1933| , Faylor 1936a , b) show clearly the existence 
of a purely hydrodynamic instability. While the flow is essen- 
tially incompressible, turbulence is generated above a critical 
Reynolds number, independent of the radial profile of angular 
momentumf|. A modern review has been given by DiPrima & 



Swinney (|1985|). Most recently [Richard & Zahn| ( |1999[ ) (here- 
after RZ) have undertaken a reanalysis of Taylor's experimen- 
tal results and, for high Reynolds number flow, have interpreted 
them in terms of a turbulent viscosity (see also Sect. 

In this contribution, we adopt the view that hydrodynam- 
ically driven turbulence can sustain the viscosity in accretion 
disks. We suggest, in Sect. ^ a viscosity prescription, the /3- 
ansatz, that represents the maximum attributable to hydrody- 
namic turbulence. We show that in the limit of low mass, 
thin disks, hydrodynamic turbulence will result in the Shakura- 
Sunyaev prescription. We then discuss the implications of the 
proposed formulation for the structure and evolution of self- 
gravitating disks, noting that even for these disks, the viscosity 
prescription differs from the a-ansatz and hence removes a dif- 



ficulty first noted by Paczyriski (1978). Finally, we discuss pro 



tostellar, galactic and galactic center disks as examples where 
the /3-ansatz may be relevant. 

2. Prescription for Turbulent Viscosity 

2.1. Reynolds viscosity as the general case 

As noted in Sect. I the need for some kind of turbulent viscos- 
ity in accretion disks is generally recognized, as is the very high 



Reynolds number of the flow in the absence of such a viscosity. 
Here, we wish to investigate in particular the case of accretion 
disks where the magnetic fields do not play an important role. 
In these circumstances, it seems reasonable to assume that the 
turbulence is driven by the velocity field in the disk, which it- 
self has characteristic length and velocity scales s (the radius 
of the orbit) and (the azimuthal velocity), respectively. 
As has been pointed out, for example, by Lynden-Bell & 



Pringle (1974) and Thompson et al. (1977), the high corre 



sponding Reynolds number 3? = sv^ / v should lead to the gen- 
eration of turbulence and hence to a steady enhancement in the 
effective viscosity. This will continue until the Reynolds num- 
ber has been reduced to approximately its critical value 5ftcrit- 
Typical values for Kdit in laboratory flows are of the order of 
^ 10^ — 10^. This limiting Reynolds viscosity can, in this case, 
be as high as 



where /3 is a constant satisfying 



1 



10"^ - 10"^ 



(2) 



(3) 



One of the cases investigated by Wendt indeed has a rotation law 
which approximates closely the profile in a Keplerian disk. 



In terms of previously introduced quantities we may write vt ^ 
(3iv^ and k (32S so that P ~ /3iP2 ~ lO"^ - "^. 

In support of this choice of s as the natural length scale 
we note that it is the only length scale which is relevant for 
angular momentum transport and which contains information 
about the driving agent for the turbulence — namely the rotation 
field; likewise, the orbital velocity is the only velocity scale 
containing such information. 

This appr oach receives furt her support from the reanalysis 
by RZ of the Wendt (1933) and Fayloi ( 1936a , b) experiments 
on turbulent viscosity generated in the flow between coaxial ro- 
tating cylinders. We note, however, that it is difficult to make 
precise comparisons between accretion disk and rotating cylin- 
ders in view of quite different constraints on the fluid flow. 

Using a definition of 3? ~ RAilAR/v appropriate to the 
experimental situation (here R is the average cylinder radius 
and Aft and AR are the relative angular velocity and gap size 
between the cylinders), RZ derive expressions for SRcrit as a 
function of relative gap size AR/R. For small gap size they 
find SfJcrit < 2000 independent of gap size. For large relative 
gap size they find that Kcrit = 5igrad(Ai?/i?)^ where their 
gradient Reynolds number ^g^iid ^ r^{dn/dR) ~ 10^ is es- 
sentially constant. Thus for small gaps the experimental data 
yield essentially the same value of SRdit as we have adopted 
in Eq. (^. For large gap sizes, the constancy of 3?grad leads to 
essentially the same functional form as in Eq. (^ but with a sig- 
nificantly smaller value of /3. RZ aiTive at similar conclusions 
for the two regimes of gap sizes from an analysis of the torques 
exerted on the cylinders in cases where the flow is turbulent. We 
believe that the small gap limit is the more relevant to the ac- 
cretion disk case, for which the speed of rotation is constrained 
at each radius by the gravitational field. The flux of angular 
momentum at each radius is thus determined by the imposed 
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orbital velocity field. By contrast, the experimental configura- 
tion constrains the flow velocity only at the inner and outer 
radius with the flow in the gap region able to take up a velocity 
profile determined by viscosity and in which the angular mo- 
mentum flux is independent of radius. In our opinion the exper- 
imental results provide strong support for a turbulent viscosity 
generated by hydrodynamic ally driven turbulence. While we 
recognize that more work is required on the question, we also 
believe that the small gap results provide significant support for 
the viscosity prescription given in Eqs. (0) and (^. 

In the following we will refer to Eq. as the /3-ansatz and to 
the disk structure arising from this viscosity prescription as f3- 
disks. We suggest that the /?-ansatz is the most appropriate ini- 
tial formulation for accretion disks since it is directly connected 
to the driving mechanism. It establishes the maximum value of 
the viscosity that can arise from hydrodynamically driven tur- 
bulence. 

The actual viscosity may, however, be limited to lower val- 
ues by such phenomena as shock dissipation of turbulent en- 
ergy if the implied turbulent velocities exceed the local sound 
speed. As we show in Sect. 2.2, this yields the a-ansatz in 



conditions relevant to these non-selfgravitating accretion disks. 
However, it leads to a different prescription in shock limited 
selfgravitating disks. In the Couette-Taylor case, all velocities 
were subsonic (the flow was essentially incompressible) so that 
no additional constraints applied. This would also be the case 
in astrophysical disks in which /3i < Cg/v^. 

2.2. a-viscosity as the limiting case for shock dissipation 
limited low mass accretion disks. 

If the accretion disk is such that the local sound speed is less 
than the turbulent velocities implied by the /3-ansatz, i.e., f3i > 
Cs/v^, we may rewrite Eq. ^ as 



vtk ~ Av^As 



(4) 



where Av^ and As are the maximum representative velocity 
and length scales allowed by local conditions. Furthermore, we 
may write 



A., = ^As ^ ^As 
OS s 



(5) 



so that a restriction on Av^ implies a constraint also on As and 
vice-versa. 

If we consider turbulent elements in a smoothed out back- 
ground gas with sound speed Cg we may impose the limit that 
the turbulent velocity will approach but may not exceed Cg. 
Thus Eq. H gives 



^As 



or 



Au0 = (cs 



As 



V4, 



-Av^ 



(6) 



(7) 



with ( a quantity smaller than but of order unity. This estimate 
of As may be interpreted as the distance a hydrodynamically 
driven turbulent element can travel before loosing its identity 
due to shock dissipation. 

In a standard geometrically thin, non-selfgravitating accre- 
tion disk (i.e., a Shakura-Sunyaev, or a-disk) hydrostatic equi- 
librium in the vertical direction impHes 



h 
s 



V4, 



Using this in Eq. 0, we find for Shakura-Sunyaev disks 

As - (h 



and hence that 



(8) 



(9) 



(10) 



with a = ('^ again not too much smaller than unity. This deriva- 
tion of the Shakura-Sunyaev scaling, starting from the assump- 
tion of a Reynolds driven turbulence, depends on the disk mass 
being negligible, i.e., a vertical hydrostatic equilibrium of the 
form of Eq. || has to apply. For selfgravitating disks, Eq. |^ 
no longer applies, and thus the functional form of v will dif- 
fer from that of Eq. [l^. Note that the upper bound to As im- 
plies approximately isotropic turbulence. This is the standard 
a-ansatz but derived from considerations of rotationally gener- 
ated turbulence. 

It is worthwhile noting that this derivation of the Shakura- 
Sunyaev prescription not only yields its functional form, v oc 
hcs, but also the order of magnitude for the scaling param- 
eter, a ^ C^, where a is close to but less than unity. This 
value is consistent with values derived by comparing a-disk 
models with observations of disks, for instance in dwarf novae 



(Cannizzo, Shafter & Wheelei 1988) 



From the above, it is clear that the viscosity in accretion 
disks depends not only on the generation of hydrodynamic tur- 
bulence but also on the limitation arising from the requirement 
that the turbulence be subsonic. It also depends on whether or 
not the disk is selfgravitating. In the following Sect., we use he 
same logic to investigate the viscosity prescription in selfgrav- 
itating accretion disks, in which turbulence is limited by shock 
dissipation. 

3. Viscosity in Thin Selfgravitating Accretion Disks 

3.1. Conditions for Selfgravity in Accretion Disks 

In the following we assume that the accretion disks are geo- 
metrically thin in the vertical direction, symmetric in the az- 
imuthal direction, and stationary. We approximate the vertical 
structure by a one zone model. Then a disk model is specified 
by the central mass A/*, the radial distributions of surface den- 
sity S(s), central plane temperature Tc(s), and effective tem- 
perature Tefi(s) or the radial mass flow rate M. The relevant 
material functions are the equation of state, the opacity and the 
viscosity prescription. 
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One can estimate the importance of selfgravity by compar- 
ing the respective contributions to the local gravitational accel- 
erations in the vertical and radial directions. 

The vertical gravitational acceleration at the disk surface 
is 2-kGYi and GAUh/s^, for the selfgravitating and the purely 
Keplerian case, respectively. Selfgravitation is thus dominant 
in the vertical direction when 



Md TTS^S 1 h 



(11) 



where M^{s) is the mass enclosed in the disk within a radius s 
and is given approximately by Md ^ tts'^S. Typical numbers 
for h/s are in the range lO^'^ "^^. 

Similar considerations lead to the condition 



Md > AL 



(12) 



for selfgravitation to dominate in the radial direction. Thus, we 
can define three regimes as follows: 

- Non-selfgravitating (NSG) disks in which A/d(s) < 
(l/2)(ft,/s)A/* (i.e., the classical Shakura-Sunyaev disks) 

- Keplerian selfgravitating (KSG) disks in which selfgravity 
is significant only in the vertical direction and which satisfy 
the constraint (1/2) (/i/s) M* < Md(s) < M* 

- Fully selfgravitating (FSG) disks which satisfy ill* < 
Md(s) 

Because A/d(s) is a monotonically increasing function of s, 
all three regimes will arise in sufficiently massive, thin (h/s <^ 
1) disks. 

3.2. Selfgravitating Disks 

In this Sect., we will review the structure of selfgravitating 
(SG) accretion disks within the framework of the assumptions 
introduced above. Compared to the standard NSG models, both 
the KSG and FSG disks require modification of the equation of 
hydrostatic support in the direction perpendicular to the disk. 
Thus, while in the standard model the local vertical pressure 
gradient is balanced by the z component of the gravitational 
force due to the central object, in the SG case we have balance 
between two local forces, namely the pressure force and the 
gravitational force due to the disk's local mass. In the KSG 
case, in the radial direction centrifugal forces are still bal- 
anced by gravity from a central mass (Keplerian approxima- 
tion), while in the fully selfgravitating case we have to solve 
Poisson's equation for the rotation law in the disk. 

For an SG disk, hydrostatic equilibrium in the vertical di- 
rection yields 



P = ttGE^ 



(13) 



(Paczyiiski 1978), where P is the pressure in the central plane 



(z = 0), E is the surface mass density integrated in the z direc- 
tion, and G is the gravitational constant. 



Since details of the thermodynamics in the z direction are 
of no particular relevance to our argument, we shall assume the 
disk to be isothermal in the vertical direction. 

Integrating the equation of conservation of angular momen- 
tum gives 



M 



(14) 



with the radial mass flow rateQ M, the rotational frequency oj, 
its radial derivative uj', and a quantity ^ allowing for the integra- 
tion constant or, equivalently, for the inn er boundary condition 
For a detailed discussion of g see, e.g., puschl & Tscharnuter 



( 1991 ), [Popham & Narayanj ( |1995[ ) and ponea & Biermann 



( 1990 . For simplicity, we set the boundary condition ^ = 
in the subsequent discussion. This does not alter the essence of 
our argument, and only changes details close to the disk's in- 
ner radial boundary, since the product s'^oj increases with s. In 
fact Eq. |l^ applies in the general case (i.e., NSG and SG); for 
Keplerian disks (NSG and KSG), we may write uj' ^ —3u!/2s. 
Finally, we have for the sound velocity 



(15) 



Eqs. |T3| and give 



2TTGm = inGph 



(16) 



where p = E/2/i is a vertically averaged mass density. 

On the other hand, the Jeans condition for fragmentation in 
the disk into condensations of radius R is 



^qGpR' > cl 



(17) 



(see Mestel 1965 ) where q is factor of order unity. 

Thus, a selfgravitating disk is on the verge of fragment- 
ing into condensations of radius R ^ h unless these are de- 
stroyed by she ar motion associated with the Keplerian veloc- 
ity field. Thus Paczyiiski ( 1978 ) and l ater K ozlowski, Wiita & 
Paczyriski ( 1979 ) and Lin & Pringle (1987) proposed that the 
viscosity prescription was directly coupled to the above gravi- 
tational stability criterion. 

To solve for the dynamic and thermal structure of SG disks, 
a viscosity prescription has to be specified. As in the NSG case 
it is possible but not necessary that viscosity is Umited by shock 
dissipation. In the absence of such dissipation, we would, as be- 
fore, expect the /3-viscosity to apply. It is instructive, however, 
to follow the logic of Sect. 2.2 in he case of SG disks in which 
turbulence is limited by shock dissipation. 



^ We choose the convention of radial mass flow rate M and radial 
velocity Vs positive for inward motion. 
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3.3. Viscosity in shock dissipation limited 
selfgravitating accretion disks 

For a SG disk (whether Keplerian or not) Eq. ^ is no longer 
vaUd so the analysis of Eqs. ^ and [l^ no longer applies. In 
physical terms, the scale height in the disk no longer reflects 
global properties of the disk (mass of and distance to the central 
star) but is set by local conditions. 

For the selfgravitating case we have approximately 



For a KSG disk, this in turn yields 



(18) 



where conditions A/» ^ M^, and ^ M^. distinguish be- 
tween the Keplerian and the FSG cases, respectively, with the 
KSG disks as an intermediate case. 
From Eq. ^ 



SCs / 2TThTtS^ 

As C— ~ C 



1/2 



C 



1/2 



(19) 



where is defined by A/* = tts^E*. 

At the transition to the NSG regime, Eqs. |ll]and |l9|and the 
condition A/, give as before 



As - Ch 



(20) 



and hence a smooth transition to the a-ansatz. 

For the FSG regime, Eq. |l9| and the condition Md 3> A/* 
give the simple asymptotic form 



As - C(2/is)'/' 

and hence a viscosity of the form 



(21) 



(22) 



where 7 is a factor of order unity. 

The situation is more complex for intermediate values of 
Md/M^,. The derived viscosity differs in all SG cases from the 
standard a-ansatz, but approaches that form as A/j/A/* — > 0. 
Thus when hydrodynamically induced turbulence is limited by 
shock dissipation, the resultant viscosity reflects local condi- 
tions and takes the standard Shakura-Sunyaev form only when 
the disk mass is negligible. We show below that this new pre- 



scription removes a problem previously noted by Paczynski 



(1978) and others, with the structure of KSG disks with a- 



viscosity. 

3.4. Structure ofSG Disks with Shock Limited Viscosity 
It foflows from Eqs. Hj |l4| and |l5[ with ^ = 0, that 

GhM / d In s ■ 



2u \d\nu}^ 
If one adopts the standard a-prescription, this yields 

GM /dins' 



2a Vdln 



(24) 



/ • \ 2/3 
2 _ / 2GM\ _ kT^ 

^ \ 3a I inn 



or 



T. = 2.42K i- 



(25) 



2/3 



a lO-SMo/yr 



(26) 



A similar result arises for the FSG case, albeit with a different 
numerical factor resulting from the solution of Poisson's equa- 
tion. 

Thus, for a SG disk, the a-ansatz leads to the requirement 
of a constant temperature for all radii s (or, if ^ 7^ in Eq. |l^ 
the temperature is prescribed as a function of s), independent of 
thermodynamics. While the exact constancy of the temperature 
may very well be an artefact of our simplified one-zone approx- 
imation for the vertical structure, there is no reason to expect 
that proper vertical integration of the structure will change this 
fundamentally. 

The a ansatz for a SG disk also requires that the disk struc- 
ture satisfy 



M2/3 



27rG 



(3a)'/^ {2ttG) 



1/3- 



(27) 



In a standard NSG accretion disk the temperature is a free pa- 
rameter which is determined by the energy released by the in- 
ward flow of the disk gas (M), by the local viscosity, and by the 
respective relevant cooling mechanisms. The viscosity depends 
on Tc via Eq. ([l|) and on the equation of hydrostatic support in 
the direction normal to the disk (Eq. (H), which in the non-SG 
case replaces Eq. ([l3|)). 

In the SG case, it is the surface density E (and hence h) 
which must adjust in order to radiate the energy deposited by 
viscous dissipation and provided by the inward flowing mate- 
rial. While detailed solutions are beyond the scope of this pa- 
per, they clearly exist formally. On the other hand, the normal 
thermostat mechanism does not operate, at least in the steady 
state. Indeed in certain circumstances, the condition of constant 
mid-plane temperature appears to be inconsistent with the basic 
thermodynamic requirement that the average gas temperature 
in the disk exceed that of the black body temperature required 
to radiate away the energy dissipated by viscous stresses (see 
Appendix). It is therefore doubtful whether a physically plau- 
sible and stable quasi-steady state solution exists. 

On the other hand, if one adopts the alternative prescription 



(23) for shock limited viscosity proposed in Sect. 3.3, the above 



problem with constant or prescribed mid-plane temperature 
disappears, the temperature once again depends on h, and the 
normal thermostat can operate. While this does not prove the 
validity of the shock limited viscosity prescription given in 
Sect. 3.3, it is certainly an interesting consequence. 
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4. Selfgravitating /3-Disks 

4.1. General Observations 

A general analysis of SG disks is complex and beyond the 
scope of this paper. In this Sect., we examine the structure of 
/3-disks, in which the turbulence is subsonic at all radii. Before 
doing so, we make the following general observations. 

First, with the /3-viscosity prescription, Eqs. |l4|and ^give 



and 



Mu 
2tivslo' 



GhM OJ 



M 



GhM 



2v suj' 2f3s^uj' 



(28) 



(29) 



Thus the SG /3-disks recover the thermostat property of the 
standard disk, namely that the temperature and scale height can 
adjust to accomodate (radiate away) the energy input to the sys- 
tem from viscous dissipation and inward motion. 

Second, we note that if the disk matter distribution is 
clumpy (e.g., clouds within a low density smoothed out dis- 
tribution) then there may be a formal connection between the 
a- and /3-prescriptions. Since in the /3-formulation the clump 
velocities are of order shock heating will tend to heat the 
low density inter-clump gas until its sound speed Cg ^ v^. The 
inter-clump gas then has a scale height h ^ s, the scale of the 
clumpy disk, and will hence be roughly a spherical structure. 
At this point the a- and /3-prescriptions look formally identi- 
cal but the scale height and sound speed now refer to a more 
or less spherical background distribution of hot gas in which a 
disk structure of cloudy clumps is imbedded. 

4.2. Keplerian Selfgravitating (3-Disks (KSG) 

For the particular case of a KSG /?-disk we have from Eq. ^ 

c2 3GM 1 



3(3 {GM^y/^s^/^ 



For the SG /3-disk it follows immediately from Eq. ( |2^ ) and 
from mass conservation in the disk that the radial inflow veloc- 
ity Vs is given by 



M 



2, / 



-(3s^uj 



2ttsY. 

For the KSG /3-disk, we then have 
3 3 



(31) 



(32) 



Thus at each radius the inward velocity is the same fraction of 
the local orbital velocity. From Eq. ( pl| ) this, in fact, holds for 
any SG /3-disk in which the angular velocity is a power law 
function of s with adjustment only to the numerical factor in 
Eq. (p2[). If (3 satisfies the constraint (^, then the approxima- 
tion of centrifugal balance in the radial direction remains well 
justified. 



Under these conditions the dissipation per unit area of a SG 
/3-disk is given by 



D 




(33) 



where a is the Stefan-Boltzmann constant. For an optically 
thick KSG disk this yields the same radial dependence of Teff 
as for the standard disk, namely 



Tcff — 



/ GMM^ 



1/4 



,-3/4 



(34) 



This temperature dependence which is identical to that of the 
standard model then leads to the well known energy distribu- 
tion for an optically thick standard disk of F^, cx v^^^. This also 
implies that — as long as the disks are not fully selfgravitating — 
it is hard to distinguish between an a- and a /3-disk model ob- 
servationally. 

4.3. Fully Selfgravitating (3-Disks (FSG) 

We turn now to the case of the fully selfgravitating (FSG) (i- 
disk, in which the disk mass is sufficiently great that it dom- 
inates the gravitational terms in the hydrostatic support equa- 
tion in both the radial and vertical directions. While there are 
many potential solutions for the FSG disk structure, one is well 
known in both mathematical and observational terms, namely 
the constant velocity (u^ = const.) disk. Within such a disk 
structure we have simultaneous solutions to the equation of ra- 
dial hydrostatic equilibrium and Poisson's equation of the form 



(35) 



Vd, = SUJ ~ vq and E oc Sq — 



( fToomrel |l963| , [Mestej |l963| ). For the FSG disk, Eq. ^ then 
leads to 



(30) E 



M 



2-k(3svq 



(36) 



which has the same radial dependence as the structural solution 
shown in Eq. (35 ). Thus Eq. (56|) may be viewed as giving the 



rate of mass flow through the disk for a FSG /3-disk with con- 
stant rotational velocity uq. Finally, the equation of continuity 
provides a constraint if the structure is to maintain a basically 
steady state structure. For the = = const, disk, this yields 



as 



d d 
— (sWsE) = T^(/3i;0sE) = 
OS OS 



(37) 



Thus the constant velocity disk represents a steady state solu- 
tion in regions sufficiently far from the inner and outer bound- 
aries of the /3-disk. 

It is then possible, in the spirit of the discussion of Eqs. 
and (^4|), to calculate the energy dissipation rate per unit area 
D for the constant velocity /3-disk. We then find 



(38) 
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so that 



off 



Mvl 



1/4 



-1/2 



(39) 



The flux density, F^, emitted by an optically thick, constant 
velocity /3-disk it then given by 



(40) 



In reality, a sufficiently massive disk may be expected to have 
an inner Keplerian (standard) zone, a Keplerian selfgravitating 
zone (KSG), and a fully selfgravitating zone (FSG). We should 
therefore expect a smooth transition in the spectral energy dis- 
tribution from the v^^'^ spectrum of the inner two zones to 
the v^^ spectrum arising at longer wavelengths from the FSG 
zone. The transition frequency z^trans may be derived by solv- 
ing Eqs. (^4|) and (^9|) for s, determining a value of transition 
temperature Ttians and setting 




(41) 



One could turn this argument around and argue that, if no other 
components contribute to the spectrum, the flatness of the vFi, 
distribution is a measure for the importance of selfgravity and 
thus for the relative mass of the accretion disk as compared to 
the central accreting object. This, of course, applies only to the 
optically thick case which may not arise frequently in strongly 
clumped disks. 

4.4. Time Scales 

The evolution of accretion disks can be described by a set of 
time scales. For our purposes, the dynamical and the viscous 
time scale are of particular interest. 

The dynamical time scale rdyn is given by 



_ 1 



(42) 



While this formulation applies to all cases, selfgravitating or 
not, it is only in the non-SG and in the KSG cases that u) is 
given by the mass of the central accretor and by the radius. In 
the FSG case, lo is determined by solving Poisson's equation. 
The time scale of viscous evolution Tvisc is given by 



(43) 



In the standard non-SG and geometrically thin {h <^ s) case (a 
disks), this leads to 



noil— SG 
vise 



\h) a 



'^dyn ^ ^dyn 



(44) 



In KSG and FSG disks (/3-disks), 

'^dyn 



is given by 



KSG 



FSG 



SG 



(45) 



With a < 1 and /3 ^ 1 (Eq. ^ under all circumstances Tvisc ^ 
Tdyn- In the SG cases the ratio between the two time scales 
decouples from the disk structure. In all cases the models are 
self-consistent in assuming basic hydrostatic equilibrium in the 
vertical direction. 

5. Possible Applications 

5.7. Protoplanetary Accretion Disks 

T Tauri stars have infrared spectral energy distributions vF^ 
which can be approximated in many cases by power laws 
vFi, oc v"" with a spectral index n in the range ^ ... 1.3. 
Assuming this spectral behaviour to be due to radiation from 
an optically thick disk, it translates into a radial temperature 
distribution Tcff oc s^^ with n = 4 — 2/g. 

An optically thick non-selfgravitating accretion disk which 
radiates energy that is liberated through viscous dissipation, 
i.e., an active accretion disk shows a spectral distribution with 
q = 3/4 or n = 4/3. This immediately excludes optically thick 
non-selfgravitating standard accretion disks as the major con- 
tributor to T Tau spectra. 



Adams, Lada & Shu (1988 ) were the first to discuss the pos- 



sibility of a non-standard radial temperature distribution with 
q ^ 3/4. Using g as a free parameter, they find that for flat 
spectrum sources, their best fits require disk masses that are 
no longer very small compared to the masses of the accreting 
stars. They already mention the possibility that the flatness of 
the spectrum and selfgravity of the disk may be related. On the 
other hand, at that time this in direct argument was the only ev- 
idence for large disk masses. Beckwith et al. ( 1990 ) in their 
survey of circumstellar disks around young stellar objects also 
find preferentially disk spectra that are considerably flatter than 
predicted by the standard optically thick disk models. For more 
than half of their objects they derive disk masses that corre- 
spond to the KSG and FSG cases. On the other hand, Matta 
(1993) proposed that flat disk spectra are the consequence of 
dusty envelopes engulfing a star with a standard disk around it. 
Recently, Chiang & Goldreich (1997) have investigated in de- 
tail non-selfgravitating passive accretion disks, i.e., disks that 
are heated by radiation from the star and re-radiate this energy. 
Depending on the details of the flaring of the disk, this can lead 
to considerably flatter spectra than expected from active disks. 

However, in the meantime, high resolution direct observa- 
tions of protostellar disks yield independent strong evidence 



for comparatively large disk masses. Lay et al. (1994), for in- 
stance, find a lower limit for the disk masses in HL Tau — one 
of the sources in Adams, Lada & Shu's sample of flat spectrum 
T Tauri stars— of - 0.02 Mq. 

We suggest that the flatness of the spectrum actually reflects 
the mass of the disk, i.e., the importance of selfgravity. For 
disk masses considerably smaller than ~ 1/30A/*, the stan- 
dard accretion disk models apply. For disks whose masses are 
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larger but still small compared to A/* the spectral behaviour 
is not altered significantly, but the disk structure and the time 
scale of disk evolution (Tvisc, see Eqs. ^ and ^5|) change. For 
even more massive disks, we expect a clear trend towards flat- 
ter spectra that approach an almost constant i^F^ distribution if 
selfgravity in the disk becomes important in the radial as well 
as in the vertical direction. 

5.2. Galactic Disks 



The relevance of viscosity in the evolution of galactic disks 
has been the subject of discussion since von Weizsacker (1943, 
1951) and j^iisl (1952) first raised the issue nearly fifty years 
ago. They noted then that, with an eddy viscosity formulation 
(a /3-disk), the time scale for evolution of typical galactic disks 
was comparable to the age of the universe and suggested that 
this might account for the difference between spiral and ellip- 
tical galaxies. 

With the subsequent realization that galactic disks moved 
primarily under the influence of extended massive halos, inter- 
est in FSG disks waned. However, as noted above, it is possible 
for a massive disk to exist and evolve under the influence of 
viscosity, while embedded in such a halo gravitational field. 
Indeed, in the event that such a structure forms, it must evolve 
under viscous dissipation and can achieve a quasi-steady state 
with essentially the same mass and energy dissipation distribu- 
tion as for the FSG constant velocity disk. We refer to this case 
as an Embedded Self-Gravitating (ESG) disk. 

The time scale for viscous evolution TvIsc as given in Sect. 
4.4 suggests a means of differentiating between the a- and 
/^-formulations for this case. For a normal spiral galaxy with 
a suggested mean temperature in the gaseous disk of around 
10"* K and a scale height of around 300 pc, we obtain 



at least slows down) in the constant angular velocity inner re- 
gions ([eke 1979, Dather & Biermann 1990). Similarly, sev- 



lOVdyn - 3 10" yr 



r^^') ~ 10^ 

Vise 



10Vdy„~310^-310i°yr 



(46) 



Thus, with these parameters, little evolution would take place 
in a Hubble time on the a-hypothesis but significant evolu- 
tion is predicted on the /3-hypothesis. This problem of the 
viscous time scale in a selfgravitating a accretion disk was 
also noted by Shlosman & Begelman (1987, 1989). Shlosman, 
Frank & Begelman (1989) proposed non-axisymmetric distur- 



bances ("bars within bars") as an alternative way of transport- 
ing angular momentum in the radial direction within a suffi- 
ciently short time scale. 

In terms of inflow velocities the /3-ansatz suggests values 
in the range 0.3 — 3 kms^^ which would be exceedingly hard 
to measure directly. The a-ansatz suggests still lower values. 
On the other hand, it may be possible to provide limits on 
the viscosity through other observational constraints. For ex- 
ample, the build up of the 3 kpc molecular ring in our own 
galaxy can be interpreted as due to viscosity driven inflow in 
the constant velocity part of the galactic disk which ceases (or 



eral authors have suggested that the radial abundance gradients 
observed in our own and other disk galaxies may be due to ra- 
dial motion and diffusive mixing associated with the turbulence 
generating the eddy visc osity (|^acey & Fal| 1985, Sommer- 
Lars on & Yoskii |1990t [Koppen| |1994 Edmunds & Green- 
how 1995, Isujimoto et al. 1995). According to these authors, 
radial inflows of around 1 kms~^ at the galactic location of the 
Sun are required for optimum fits to the abundance gradient 
data within the context of the viscous disk hypothesis. Such in- 
flow velocities are consistent with the /3-ansatz but could, of 
course, be generated also by other means (e.g effects of bars, 
magnetic fields). 



5.3. Ultraluminous galaxies 

Recent high resolution imaging of ultraluminous galaxies in 
the near infrared and mm wavelengths bands shows dense gas 
and dust accretion disks in their galactic nuclei. The two nu- 
clei in the merger galaxy Arp 220, for instance, have masses 
of the order of several 10^ within radii of <^ 100 pc 
(Scoville 2000| ). Similar properties, albeit less well resolved 
than in Arp 220, seem to be typical for this class of galax- 
ies dSolomon et al.l 119971, pownes & Solomon! 119981). Most, 



if not all, ultraluminous galaxies seem to be merging galax- 
ies ( [Sanders & Mirabel 1990- In Arp 220, these gas masses 
are t he major contri butor to the dynamical mass in the two nu- 
clei ( Scoville 200C), i.e., these nuclear disks are selfgravitat- 
ing. Most likely this is true for the nuclear disks in other ultra- 
luminous galaxies as well. The merger process is presumably 
responsible for transporting large amounts of material into the 
central few hundred parsecs, thus filling a mass reservoir which 
is then available for subsequent disk accretion to the very cen- 
ter. 

Within the framework of /3-disks, one finds that the vis- 
cous accretion time scale Tvisc increases towards larger radii 
as long as the surface density E in the disk increases with 
radius s not steeper than S cx s, which is most likely ful- 
filled. Then the viscous time scale at the disk's outer edge is 
an upper limit to its evolution time scale. For Arp 220 (disk 
mass ~ 210^Mq; outer radius ^ lOOpc) one finds a time 
scale of - lO^-^ yr//3 - 3 ... 30 lO"^ years for [3 = lO^^.-.-a^ 
which, in turn yields accretion rates M ^ lO^ "^ Mq yr^^. 
Such rates lead to accretion luminosities iaccr = rjAIc^ up to 
^ T^IO^^-^- -^^-^ ergs^^, where rj 0.1) is the conversion effi- 
ciency of gravitational energy into radiation and c is the speed 
of light. Such luminosities are large enough to power even the 
strongest AGN and the time scales are very much shorter than 
the Hubble time. 

Assuming that these rates can be maintained during a size- 
able fraction of Tvisc, a significant fraction of the disk's original 
gas mass could be accreted to much smaller radii, presumably 
to a black hole in the very center (some will be lost to star for- 
mation or winds). In this process the black hole gains a consid- 
erable amount of mass within a relatively short time scale. One 
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may speculate that this is actually the process that produces the 
most massive black holes in the young universe. By contrast, 
galaxies that do not undergo mergers presumably have no way 
of rapidly collecting such large masses of gas within 10^ pc. 
As a consequence, these disks are less likely to be selfgravi- 
tating and thus are likely to have longer Tvisc- The nuclei of 
such galaxies will accrete much smaller amounts of material 
over longer time scale, resulting in lower mass central black 
holes ( puschl 1988a, b). An example may be our own Galactic 
Center 

6. Summary 

We propose a viscosity prescription based on the assumption 
that the effective Reynolds number of the turbulence does not 
fall below the critical Reynolds number In this parametriza- 
tion the viscosity is proportional to the azimuthal velocity and 
the radius (/3-disks). This prescription yields physically consis- 
tent models of both Keplerian and fully selfgravitating accre- 
tion disks. Moreover, for the case of thin disks with sufficiently 
small mass, we recover the a-disk solution as a limiting case. 

Such /3-disk models may be relevant to protoplanetary ac- 
cretion disks as well as to galactic and galactic center disks. 
In the case of protoplanetary disks they yield spectra that are 
considerably flatter than those due to non-selfgravitating disks, 
in better agreement with observed spectra of these objects. In 
galactic disks, they result in viscous evolution on time scales 
shorter than the Hubble time and thus offer a natural expla- 
nation for an inward flow that could account for the observed 
chemical abundance gradients. In galactic centers, /3-disks may 
be the supply for powering AGN and for forming supermassive 
black holes within time scales short compared to the Hubble 
time. 

Finally, /3-disks yield a natural solution to an inconsistency 
in the a-disk models if the disk's mass is large enough for 
selfgravity to play a role. This problem arises even in Keple- 
rian selfgravitating disks in which only the vertical structure is 
dominated by selfgravity while the azimuthal motion remains 
Keplerian. 



= 8.53 10^ K 



1/4 ^1/4^-3/4 
1/4 



(A.2) 



with m the mass of the central star in solar units and sa the 
radius in astronomical units. 

An essential thermodynamics requirement is that T^. > T^g 
or that 



1^-3.5310- ^ 



a_!^ < 1 



(A.3) 



This condition is satisfied provided that 



m > mx = 3.2710 



-4' 



9/5 



(A.4) 



and that the disk is selfgravitating in the vertical direction at 
SA. The latter condition leads to a second requirement on m. 

For a standard Keplerian disk, the mass flow rate is given 
(Eqs.limpby 



M « — « -4r- = aMd {-] oj (A.5) 



with A/d the disk's mass. From Eq. [rT|, the condition that the 
disk is non-selfgravitating is il/d < {h/2s)M^, and hence, from 
Eq. O, that 



M < 



3 /^A.r3\ 1/2 



a fhy fGM: 



2 \ s J \ s3 



(A.6) 
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Appendix A: Thermodynamic Considerations for KSG 
a-disks 

For a KSG a-disk, we have from Eq. ^ that 

\ 2/3 



T, = 2.42 K 



1 



M 



a 10-6Mo/yr 



(A.l) 



If the disk is optically thick and advection is negligible, viscous 
dissipation leads to local effective temperature of 



m < rfiQ ~ 3.14 10 



a_i 1 - 



»3/2 
3/2 



(A.7) 



A selfconsistent and physically acceptable solution can be ob- 
tained only if niQ > tot, that is the disk becomes selfgravitat- 
ing at values of to which are sufficiently high that thermody- 
namic requirements are not violated. This condition may then 
be written as 



> 1.01 10" 



-3/10 -1/10 1/5 



(A.8) 



Thus thin KSG a-disks appear to be inconsistent with basic 
thermodynamic requirements if m 1, sa ^ 1- There is no in- 
consistency if either or both of these quantities are sufficiently 
large. 
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